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We present the significant detection of the first extragalactic pulsar wind nebula (PWN) detected in gamma rays, N 157B, located in the large 
Magellanic Cloud (LMC). Pulsars with high spin-down luminosity are found to power energised nebulae that emit gamma rays up to energies 
of several tens of TeV. N 157B is associated with PSR J0537-6910, which is the pulsar with the highest known spin-down luminosity. The High 
Energy Stereoscopic System telescope array observed this nebula on a yearly basis from 2004 to 2009 with a dead-time corrected exposure of 
46 h. The gamma-ray spectrum between 600 GeV and 12 TeV is well-described by a pure power-law with a photon index of 2.8 ± 0.2s,a, ± 0.3sys, 
and a normalisation at 1 TeV of (8.2 ± 0.8siai ± 2.5syst) x 10"" cm^'s^'TeV"'. A leptonic multi-wavelength model shows that an energy of about 
4 X 10*" erg is stored in electrons and positrons. The apparent efficiency, which is the ratio of the TeV gamma-ray luminosity to the pulsar's spin- 
down luminosity, 0.08% ± 0.01%, is comparable to those of PWNe found in the Milky Way. The detection of a PWN at such a large distance is 
possible due to the pulsar's favourable spin-down luminosity and a bright infrared photon-field serving as an inverse-Compton-scattering target 
for accelerated leptons. By applying a calorimetric technique to these observations, the pulsar's birth period is estimated to be shorter than 10 ms. 

Key words, gamma rays: general - pulsars: individual: PSR J0537-6910 - supernova remnants: individual: N 157B - Magellanic Clouds 



1. Introduction 

In recent years, many Galactic pulsar wind nebulae (PWNe) 
have been dis covered to be gamma-r ay emitters (for a re- 
Send offprint requests to: Nukri Komin (komin@lapp.in2p3.fr) view, see e.g. Ide Ona Wilhelmil |201 ih . It was predicted by 
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lAharonian et alJ (Il997h that the gamma-ray luminosity of these 
nebulae is connected to the spin-down power E, i.e. the 
loss rate of rotational energy of the pulsar, and that pul- 
sars with E > 10-'"^""'^((i/lkpc)^ erg s"' (with d being the 
distance to the pulsar) power nebulae that are detectable in 
gamma rays. On the basis of the Galactic Plane Survey carried 
out by H .E.S.S. (|Aharonian et al. 2005; Aharo nian et al. 2006b: 
[past et al. 2012 h , Carrigan et al. (2008) suggested that pulsars 
with a spin-down power of £ > 10"^'*(t//lkpc)^ erg s"' may be 
correlated with nebulae detectable by H.E.S.S. 

The most energetic pulsar known is PSR J05 37-6910, with 
a spin-down power of £ - 4.9 x 10^^ erg s ' (M arshall et alj 
Il998t [Mancheste r et alJl20051 PSR J0537-6910 is also one the 
most distant pulsars known to date: located in the LMC, it has 
an estimated distance of (48.1 ± 2.3stat ± 2.9syst)kpc dMacri et"al] 
[2OO6). This pulsar is the compact central object of the super- 
nova remnant (SNR) N 157B (also called LHA 120-N 157B, 
SNR B0538-691, or NGC 2060).N157B has been observ ed ex- 



tensively in X-ra ys (IWang & Gottheljll998l: [Wang et a l] 1200 it 
Chen et alj20 0^. N 157B is a Crab-like SNR fcott helf & Wand 



19961) 7 the emission from N 157B is dominated by synchrotron 
radiation from the PWN, which shows a bar-like feature sur- 
rounding the pulsar representing the reverse shock of a toroidal 
wind from the pulsar as well as a long tail of diffuse emis- 
sion of about 20" X 30" in the north-west direction from the 
ba^ (Wang et al. 2001; Chen et al. 2006) (see also Fig. Wp). 
IChen et alj ( l2006l) show that faint thermal emission with a diam- 
eter of about 100" is observed from the supernova ejecta and, un- 
like shell-type or composite SNRs, no emission from the super- 
nova forward shock is detected. The supernova ejecta apparently 
expand into the low-density interi or of the superbu bble formed 
by the stellar association LH99 dChen et all 120061). The char- 



acteristic age of the pulsar is about 5000 years ( Marshall et alj 
I 998h. which is consis tent with estimates of the SNR's age 



( Wang&Gott"hell[T998 ). ;^arge-scale diffuse emission from the 
gen eral direction of N 157B was discovered with the Fermi satel- 
lite (lAbdo et al.ll2010l) . which was interpreted as emission from 
massive star-forming regions. No significant emission, in addi- 
tion to the diffuse emission, has been detected from N 157B. 
Despite its extreme distance, PSR J0537-6910 has a large 



enough spin-down flux of E/d^ = 
a nebula detectable with H.E.S.S. 



2x10^^^ ergs 'kpc 



to power 



minimisation. The remaining background was estimated from 
both rings around each sky position to generate the gamma-ray 
image and spatial analysis (ring background, Berge et al. 200_3) 
and test regions with similar offsets from t he camera centre f or 
the spectral analysis (reflected background, iBerge et al]|2007h . 

Figure[T^ shows a gamma-ray excess image of N 157B. The 
image was smoothed with the H.E.S.S. point-spread function, 
where 68% of the events were contained in a circle with a radius 
of 0.06°. The significance was calculated from the counts in a 
circular integration region with a radius of 0.06° around each 
sky bin. Fitting a point-like source folded with the instrument's 
point spread function results in a best-fit position of the source 
of RA = 5''37™44^ Dec = -69°9'57", equinox J2000, with a sta- 
tistical uncertainty of + 1 1" in each direction; the source is hence 
labelled HESS J0537-691. With the H.E.S.S. standard pointing 
correction, point-like sourcesj;an be localised with a systematic 
uncertainty of 20" per axis dAcero et al. 20101) . FigurefT}? shows 
an X-ray image of the supernova remnant and its central PWN 
with overlaid confidence contours of the gamma-ray source po- 
sition. The best-fit position is consistent with the pulsar position; 
the slight offset from the pulsar along the tail of the PWN is not 
significant (1.3 cr of the combined statistical and systematic er- 
ror). 

In a circular region with a radius of 0.1° around the pul- 
sar position, 395 gamma-ray candidate events were found. A 
total number of 3152 events were found in a background re- 
gion with an area that is larger by a factor of 18.62. The cor- 
responding gamma-ray excess is 226 events with a statistical 
sign ificance of 14 cr (calculated using formula (17) of Li&]\^ 
11983 ). The photon spectrum between 600 GeV and 12TeV of 
the gamma-ray excess can be described by a pure power-law, 
dN/dE = ^o(E/lTeV)'^, with a normaHsation at 1 TeV of 

00 = (8.2 + 0.8.stat) X lO^'^cm^^g-i'TeV-' and a photon index 
of r = 2.8 + 0.2stat, which corresponds to an energy flux between 

1 TeV and lOTeV of (1.4 ± 0.1) x lO'^^ergcm' ^ s"' or 2% of 
the energy flux of the Crab Nebula ( lAharonian et al. 2006 a). The 
gamma-ray luminosity between 1 TeV and 10 TeV for a distance 
of48<i48kpc is (3.9 + 0.3) x IQ^^dj^ergs-^ This resuh was de- 
rived from data accumulated over 4 years during the summer 
rainy seasons with varying atmospheric conditions. From the 
analysis of different subsets of the data, a systematic uncertainty 
of 30% in the flux and 0.3 in the photon index is estimatecfl. 



2. High Energy Stereoscopic System observations 
and results 

The High Energy Stereoscopic System (H.E.S.S.) 
(lAharonian et al. 2 006a) is a system of four Imaging Cherenkov 
Telescopes, located in the Khomas Highland of Namibia at an 
altitude of 1800m. The H.E.S.S. telescope array's location in 
the southern hemisphere is ideal for observing the Magellanic 
Clouds. Furthermore, H.E.S.S. is sensitive to gamma rays at 
energies above 100 GeV up to several 10 TeV. The arrival 
direction of individual gamma rays can be reconstructed with 
an angular resolution of higher than 0.1°, and their energy is 
estimated with an relative uncertainty of 15%. 

The region containing N 157B was observed with H.E.S.S., 
the data presented having been recorded from 2004 to 2009, with 
a dead-time corrected exposure of 46 hours after data quality se- 
lection. The observations were carried out at large zenith angles 
(with a mean of 47°) leading to an elevated energy threshold 
of about 600 GeV. The data wer e analysed using Model analysis 
with standard cuts dde Naurois & Rolland 2009), where the cam- 
era images are compared with simulations using a log-likelihood 



3. Discussion 

The observed gamma-ray luminosity between 1 TeV and 10 TeV 
of the PWN con-esponds to (0.08 + O.Ol)dl^% of the spin-down 
power o f the pulsar, a value typic al of young PWNe (for a review 
see e.g. Ide Ona Wilhelmill20ll . F^ilsar wind nebulae are very 
efficient in producing TeV gamma rays via the inverse Compton 
(IC) up-scattering of lower-energy photons by relativistic elec- 
tronfl This makes them not only the most abundant source 
type among the TeV emitters in the Milky Way, but also with 
HESS J0537-691 the first extragalactic TeV source that is un- 
related to a galaxy or an active galactic nucleus. N 157B can be 
detected at such a large distance not only due to the high spin- 
down power of the pulsar but the strong infrared photon-fields 
from nearby sources serving as additional targets for the IC scat- 
tering. 



' These uncertainties are specific to this data set and should not be 
used for any other H.E.S.S. result. 
- The term electrons refers to electrons and positrons. 
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Fig. 1. Gamma-ray and X-ray view of N 157B. a) Image of the TeV gamma-ray excess obtained with H.E.S.S. The image was 
smoothed with the instrument's point-spread function, which is shown in the inset. The green contour lines denote 6, 9, and 12 cr 
statistical significances. The black box outlines the region shown in the X-ray image, b) Image of the X-ray flux in the energy band 
from 0.8 keV to 8 keV {Chandra ACIS-S, observation n°2783, 48.2 ks dead-time corrected exposure). The data were smoothed with 
a Gaussian of cr = 1 .5" and the exposure map was computed with an assumed spectral index of 2.3, which is the average index of the 
nebula. The colour coding represents the X-ray flux in photons cm"^ s ; for better visibility, it is clipped at 10"'^ photons cm"^ s"' 
at the lower end and saturates at 10"^ photons cm"^ s"'. The white contour lines denote regions of 68%, 95%, and 99% confidence 
for the position of the gamma-ray source HESS J0537-691. The position of the pulsar PSR J0537-6910 is indicated by a black 
cross. The green cross marks the best-fit position and the systematic uncertainty in the pointing of the H.E.S.S. telescopes. 



Figure |2] shows the spectral energy distribution (SED) of 
N 157B. It exhibits synchr otron emission that ha s been detected 
at radio wavelengths (see 'Lazendic et al."200C i, and referen ces 
therein) and in X-rays (Chen et al. 2006; Denne rfet al.ll200ll) . as 
well as TeV gamma-ray emission from IC scattering. In addition 
to the cosmic microwave background (CMB), far infrared pho- 
tons from the OB association LH 99 and the nearby star-forming 
region 30 Doradus are imp ortant targets for IC scat tering. Using 
observations from Spitzer (llndebetouw et al.ll2009l) . the infrared 
photon fields are modelled as black-body radiation with a tem- 
perature of 80 K and an energy density of 8 .9 eV cm""* for LH 99, 
and a temperature of 88 K and an energy density of 2.7 eV cm"^ 
for 30 Doradus. These are only upper limits to the infrared fields, 
as the (unprojected) distances of N 157B to these objects are un- 
known. 

There is no evidence of extended emission beyond the angu- 
lar resolution of H.E.S.S., thus the synchrotron and IC emission 
regions cannot be separated spatially. Hence, a simple one-zone 
model is assumed, where only a single electron population is re- 
sponsible for both the synchrotron and IC emission. An electron 
spectrum following a broken power-law with two breaks and 
covering the energy band from 1 eV to 10'^ eV is adopted. The 
low energy b reak is an intrinsic break of the injection spectrum 
proposed by IVenter & de Jagerl (l2007h : the high energy break 
arises from the cooling of the particles. It is assumed that the 
cooling break appears at the energy where the synchrotron loss 
time (which depends on the magnetic field) is equal to the age of 
the remnant (which is assumed to be 5000 years). A low-energy 
spectral index of 1.33 is used to reproduce the radio spectrum, 
and a spectral index of 3.60 above the cooling break is consis- 
tent with both the X-ray and TeV data. A ssuming that the cool- 
ing steepens the spectral index by one (IKardashevI 1 1 962h . the 



uncooled high-energy spectral index is 2.60. This model, which 
is represented by the green lines in Figure |2] requires an intrin- 
sic break at 54GeV and a magnetic field of 41 yuG, the coiTe- 
sponding cooling break being at 1 . 1 TeV. Assuming a distance 
of 48 kpc, the total energy Wtot stored in electrons is 4 x 10'*'' erg. 
In a second, more conservative model, any assumptions about 
the age and cooling of the particles are abandoned to minimise 
the total energy stored in electrons. The intrinsic break is set to 
7 GeV, the lowest energy still being compatible with the radio 
data, and the cooling break is set to 4.22 TeV, the highest energy 
being compatible with the TeV data points. These prerequisites 
require an uncooled high-energy spectral index of 2.57 and a 
low-energy spectral index of 0.93 to reproduce the radio data. 
This model is represented by the red lines in Fig. |2j it requires a 
magnetic field of 47 /jG and the total energy is only 50% of the 
energy in the first model. The uncertainty in the distance mea- 
surement adds an error of 15% to the estimate of Wtot- Major 
uncertainties in the estimation of Wtot lie in the uncertainties in 
the gamma-ray spectrum. Varying the TeV flux by 30% while 
fixing the radio and X-ray points changes the total energy by 
25%. The total energy content Wtot cannot be determined from 
the radio and X-ray data alone, since the strength of the syn- 
chrotron emission is governed by the a priori unknown mag- 
netic field. Using the observed TeV spectrum presented here, it 
is possible to derive the magnetic field in the PWN and thus the 
energy content in electrons in the PWN. In a hadronic scenario, 
gamma-ray emission is produced in the decay of tt" mesons pro- 
duced by inelastic interactions of accelerated protons with am- 
bient material. The observed TeV emission can be described by 
the emission of a proton population following a power-law with 
an index 2 and an exponential cut-off at 23 TeV (blue, long- 
dashed line in Fig.|2]l. The total energy in this proton population 
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is Whad = 1-8 X 10^ («/cm ) 'c/^gCrg. An ambient density of 
at least 100 cm"-' would be necessary to produce this emission 
by a single supernova. That the SNR is expected to expand into 
the low-density interior of a superbubble makes this scenario un- 
likely. 

As proposed by Ide Jageii ( l2008l) . Wtot can be used to esti- 
mate the birth period of the pulsar PSR J0537-6910is relatively 
young, its characteristic age being shorter than the cooling time 
of most of the electrons. Particles have survived since the ear- 
liest epoch when the pulsar's spin was close to its birth period. 
Therefore, Wtot can be related to the pulsar's birth and current 
period by the calorimetric expression 



«7 (£'rot,0 - £'rot) 



1 

eri-I 



= 2 X 10^°e77 



10^5 gcm2 



10 ms 



10 ms 



2\ 



(1) 

(2) 

erg, (3) 



where Pq is the birth period of the p ulsar, P is the current pe- 
riod of 16.1 ms ([Marshall et al.|[T998h . rj denotes the conversion 
efficiency of spin-down power into accelerated electrons, and 
e is the relative average energy-loss rate, which takes into ac- 
count the energy already radiated by the particles during ear- 
lier epochs and adiabatic losses of the energy. For the pul- 
sar's mome nt of inertia, the canonical value of / = 10"*^ gcm^ 
is adopted (' Lattimer & Prakashll2001h : the moment of inertia 
could be higher owing to variations in the pulsar mass be- 
tween 1.4 solar masses (Mp) and 2.5 as proposed by 
iBelczvnski & TaamI j2008h Models indicate that // can be as low 
as 0.3 (ISchocketalJl20ld for MSH 15-52) or as high as 0.7 
dde Jager et al .''2008'. for G21.5-0.9). Magnetohydrodynamical 
simulations Ide Jager et al. 2009) show that the adiabatic energy 
loss is around e ~ 0.5. Further radiation losses during earlier 
epochs of the nebula reduce e to less than 0.5. Choosing 0.7 and 
0.5 for T] and e, respectively, shows that the birth period must 
have been shorter than 10ms. Using the conservative model, 
the birth period increases by 23%. This result is consistent with 
earlier estimations of the pulsar's birth period, [Marshall et alJ 
(Il998 l) estimated the birth period by comparing the pulsar's char- 
acteristic age with the age of the SNR for different braking in- 
dices, for SNR ages of more than 4000 years the birth period is 
shorter than 10 ms for a larg e range of braking ind ices. From the 
extrapolation of glitch data, [Marshall et alj (|2004|) derive a pul- 
sar period of 1 1 ms 5000 years ago. In the present paper, for the 
first time the birth period of a pulsar is obtained directly using a 
calorimetric technique, which depends on neither the glitch his- 
tory nor the braking index and is — for the conservative model 
— completely independent of the age of the remnant. 

This result confirms that PSR J0537-6910 has with a birth 
period of shorter than about 10 ms, the shortest birth period ever 
derived for a pulsar These short rotation periods are only known 
for millisecond pulsars that have been spun up after their birth 
by a companion star Simulations show that pulsar birth peri- 
ods can be related to some parameters of the progenitor stars. 
[Heger et al. ( 2005) show that more massive progenitor stars pro- 
duce heavier and more rapidly rotating pulsars: stars of 15, 20 
and 35 M© are required to produce pulsars with 1 1 ms, 7 ms 
and 3 ms, respectively. This is consistent with an earlier estimate 
of the PSR J0537-6910 progenitor mass of 20-35 Mq based on 
the comparison of the observed metal abundances i n the super- 
nova ejecta of SNR N 157B with supernova models dChen et alj 
12006) . Such massive stars are close to the threshold for the for- 
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IC emission on infrared from 30 Dor 
IC emission on CMB 
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Fig. 2. Spectral ener gy distribution of N 157B. The data shown 
are radio emissi on ( Lazendic et al.l 12000^ blue triangles), non- 
thermal X-rays dChen et al.l 20061 blue band), and TeV gamma 
rays (H.E.S.S., this work, black circles). The green, thick lines 
represent the model with a cooling break depending on the mag- 
netic field and the age of the remnant. The red, thin lines de- 
note the conservative model that minimises the energy content 
of electrons. The blue, long-dashed line represents a hadronic 
scenario. 



mation of black holes in the supernova explosion. lFryer| dl999h 
show, for instance, that 25 M p is rough l y th e limit for black hole 
formation. On the other hand. lOtt et al.l (1200 6) show that the pul- 
sar's birth period is rather unrelated to the progenitor's mass but 
roughly linearly dependent on the initial central iron-core spin. 
Birth periods of shorter than 10 ms require initial iron-core pe- 
riods of shorter than about 8 s. A very massive and/or rapidly 
spinning progenitor star therefore appears to be required to pro- 
duce a neutron star with a birth period as short as 10ms. 

Alternatively, the pulsar could be part of a binary system that 
has been spun up by its companion star, a scenario typically as- 
sumed for millisecond pulsars. A very massive and rapidly ro- 
tating star at a distance of about 0.8' from PSR J0537-6910 has 
been identified (Duftonet al. 2011). It was proposed that both 
objects were part of a binary system where mass was transferred 
from the pulsar's progenitor to its companion and that both stars 
experienced radial velocity kicks in the supernova explosion. 
Nonetheless, in this scenario, the general picture of a very mas- 
sive star producing a rapidly spinning neutron star remains un- 
changed. 



4. Conclusions 

Our principal conclusions are as follows: 

1. Gamma-ray emission from the PWN in N157B was dis- 
covered with H.E.S.S. observations. The energy flux be- 
tween ITeV and 10 TeV is (1.4 + 0.1) x 10"'^ erg cm"^ ^-i 
Located in the LMC at a distance of 48kpc, this is the 
most distant PWN ever detected in gamma rays and is the 
first individual stellar extragalactic TeV gamma-ray source. 
The PWN is powered by the most energetic pulsar known: 
PSRJ0537-6910. 

2. The TeV photon spectrum, in connection with radio and X- 
ray measurements, can be described with a one-zone leptonic 
model. From this model, the total energy stored in electrons 
in the nebula can be estimated to be 4 x lO'^^erg. 
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3. For the pulsar to provide this energy from its rotational en- 
ergy loss, the pulsar's birth period must have been shorter 
than about 10 ms. This is the shortest birth period ever in- 
ferred for a pulsar. In an alternative scenario, the pulsar might 
have been spun up by a companion star in a binary system. 

4. Assuming a direct connection between the pulsar's birth pe- 
riod and the mass of the progenitor star, the progenitor must 
have had a mass of at least 15 Mq. This is close to the 
limit for black hole formation. The pulsar PSR J0537-6910 
is therefore at the upper mass limit for neutron star produc- 
tion. 
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